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Abstract 

Using the natural connection equivalent to the SU(2) Yang-Mills instanton 
on the quaternionic Hopf fibration of S 7 over the quaternionic projective space 
HP 1 ~ S 4 with an SU(2) ~ S 3 fiber the geometry of entanglement for two 
qubits is investigated. The relationship between base and fiber i.e. the twisting 
of the bundle corresponds to the entanglement of the qubits. The measure of 
entanglement can be related to the length of the shortest geodesic with respect 
to the Mannoury-Fubini-Study metric on HP 1 between an arbitrary entangled 
state, and the separable state nearest to it. Using this result an interpretation of 
the standard Schmidt decomposition in geometric terms is given. Schmidt states 
are the nearest and furthest separable ones lying on, or the ones obtained by 
parallel transport along the geodesic passing through the entangled state. Some 
examples showing the correspondence between the anolonomy of the connection 
and entanglement via the geometric phase is shown. Connections with important 
notions like the Bures-metric, Uhlmann's connection, the hyperbolic structure for 
density matrices and anholonomic quantum computation are also pointed out. 
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I. Introduction 



Since the advent of quantum computation the importance of quantum en- 
tanglement cannot be overestimated. Maximally entangled states have made their 
debut to physics via Bell-type inequalities [1] exemplifying measurable differences 
between classical and quantum predictions. Recently entangled states have become 
important via their basic use in quantum computation processes [2], teleportation 
[3], dense coding [4], and quantum key distribution [5]. In the light of such appli- 
cations it has become evident that quantifying entanglement and understanding 
its geometry is a problem of basic importance. Efforts have been made to quan- 
tify entanglement by introducing suitable measures for it [6,7] . These approaches 
emphasize the difference between entangled and separable states by introducing 
measures usually related to the entropy of the states [8]. At the same time some 
authors have pointed out correspondences between the notion of entanglement and 
the basic geometry of the space of states . The space of states for spin-like systems 
with the composite Hilbert space H ~ C n is CP n_1 the n — 1 dimensional com- 
plex projective space. Here the slicing of the space of states for submanifolds of 
fixed entanglement was introduced and illustrated [9,10]. It was also shown how 
algebraic geometric ideas can be used to study two qubit entanglement within 
the framework of geometric quantum mechanics [11]. Recently in an interesting 
paper Mosseri and Dandoloff [12] used the quaternionic Hopf fibration to take 
another look at the problem of characterizing the geometry of two qubit entangle- 
ment. Their results have been generalized to three qubits by using the next Hopf 
fibration based on octonions [13]. 

This paper can be regarded as a further development in understanding en- 
tenglement in geometric terms. For illustrativ purposes we take the simplest two- 
qubit case and use the convenient quaternionic representation [12] of two qubit 
entanglement. In this picture the two qubit Hilbert space is fibered over the four- 
sphere S 4 which is isomorphic to HP 1 the one dimensional quaternionic projective 
space. This four-sphere is sliced to submanifolds of fixed entanglement. Our key 
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idea is the observation that the quaternionic Hopf bundle can be equipped with a 
natural connection enabling a geometric means for comparing states belonging to 
submanifolds of different entanglement. This connection provides a splitting for 
vectors corresponding to entangled states to parts representing their horizontal 
and vertical components. Using this splitting a natural metric (the Mannoury- 
Fubini-Study metric) is induced on HP 1 which is essentially the standard metric 
on the four spehere S 4 expressed in terms of stereographically projected coordi- 
nates. The geodesic distance with respect to this metric provides a natural tool 
for quantifying entanglement. This simple picture gives a further understanding of 
the results of Brody and Hughston [11] quantifying entanglement by the geodesic 
distance between the entangled state in question and the nearest separable state 
with respect to the standard Fubiny-Study metric on CP . The important new 
ingredient of our paper is the possibility of using the non-Abelian geometric phase 
(the anholonomy of the natural connection) in obtaining a further insight to the 
geometry of two qubit entanglement. Our approach beeing interesting in its own 
right also gives an interesting application of the idea of holonomic quantum com- 
putation [14]. 

The organization of this paper is as follows. In Section II. we briefly sum- 
marize some basic background material concerning two-qubit entanglement. In 
Section III. we reformulate the results of Mosseri et.al. on the Hopf fibering of the 
two qubit Hilbert space in a formalism convenient for our purposes. In Section IV. 
we introduce our geometric structures, the connection and the metric. In Section 
V. using this formalism we show that the geodesic distance between an entangled 
state \1/ and the closest separable state is a convenient quantity characterizing 
entanglement. In fact this quantity is expressed in terms of the concurrence of 
the entangled state in question. The nearest and furthest separable states ob- 
tained by parallel transport with respect to the instanton connection are just the 
ones appearing in the Schmidt decomposition of In Section VI. some examples 
showing the correspondence between the anholonomy of the natural connection , 
anholonomic quantum computation and entanglement via the geometric phase is 
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shown. In Section VII. some connections with important notions like the Bures- 
metric and Uhlmann's connection for density matrices are also given. We also 
relate these notions to the hyperbolic structure of the space of reduced density 
matrices. The conclusions and some comments are left for Section VIII. 



II. Two-qubit entanglement 

In this section we summarize well-known results concerning two-particle pure- 
state entanglement. Although formulas below are valid for wave functions of both 
particles belonging to a finite N dimensional Hilbert space H N , we have in mind 
the N = 2 case i.e. two qubits. As a starting point we write the two-particle wave 
function as 

N-l 

I*) = -Ttt E c «p\ a Ph \<*P) = ® \P) 2 , (!) 

V N a,f)=0 

where \a) 1 and \(3) 2 are orthonormal bases for subsystems 1 and 2 and = 1. 

The pure state density matrix of the total system is p = The reduced 

density matrices p\ and p2 characterizing the state of the system available to 
an observable capable of performing local manipulations merely on subsytem 1 
respectively on subsystem 2 are given by the expressions 

Pi = Tr 2 p = ^CC\ p 2 =Tr lP =±CTc, (2) 

where it is understood that p\ and p2 are N x N matrices expressed in the base 
|a) 1 and \(3) 2 respectively. The von Neumann entropy corresponding to the i-th 
subsystem is defined as 

St = -Trp4og 2Pi = - J2 X ^ )lo ^ ( n\ 2 = 1,2. (3) 

Using the fact that CC^ and C^C are hermitian matrices having the same real 
nonzero eigenvalues one can see that Si = S 2 . 

4 



For maximally entangled states we have pi = P2 = jjl where / is the N x N 
identity matrix. In this case from (2) and (3) it follows that C G U(N) and 
S = log 2 iV, hence in particular for two qubits we have C G U(2) and S = 1. For 
separable states we have C a p = X a Yp, hence in this case can be written in the 
product form \ip) <g> \ip). One can readily show that is separable if and only if 
S = 0. Since for separable states the matrix C is a dyadic product of two vectors 
the partially traced density matrices and C are all of rank one. For the two-qubit 
N = 2 case it means that |^) is separable if and only if detC = 0. Separable and 
maximally entangled states are extremal in the sense that they give the minimum 
and maximum values for the von Neumann entropy. In between these cases lie 
states of intermediate entanglement characterized by the values < S < log 2 iV. 

It is well-known that an arbitrary state G H N <8> H N expressed as in (1) 
can be transformed to the Schmidt form [15] 



by means of local unitary U(N) x U(N) transformations acting independently on 
the two subsystems. The nonnegative real numbers Xj are the eigenvalues of the 
reduced density matrices, hence they sum to one in accordance with the property 
T^Pi,2 = 1- Notice that in the sum only the nonzero eigenvalues contribute which 
are the same for both reduced density matrices. The orthonormal states 2 can 
be obtained by finding the eigenvectors corresponding to the nonzero eigenvalues 
of the reduced density matrices p\^- 

Let us give explicit expressions for the N = 2 case! We write 



N-l 




(4) 



3=0 



V2 



1 





Let us define the complex numbers 



z = ac + bd, w = ad — be, C = ab + cd. 
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(6) 



Then the reduced density matrices are 



1/H 2 + |6| 2 z \ l/|o| 2 + |c| 2 C N m 

Due to the normalization condition (^l^) = 1 we have Trpi 2 = 1, moreover 

1 2 1 2 

by virtue of (2) Detpi^ = ||DetC| = -^\w\ . The magnitude of w is the concur- 
rence satisfying the relation < C = \w\ < 1 [7]. It is obvious that for separable 
states one has C = 0. For maximally entangled states C G £7(2) ~ £7(1) x SU{2) 
hence C = |DetC| = 1. The eigenvalues of the reduced density matrices in terms 
of the concurrence read as 



x± = l(i±VT^). (8) 

Using this result it is reasonable to define the entanglement of a two-qubit pure 
state |^) to be its von Neumann entropy [16] i.e. 



= -Trp!log 2 pi = -Trp 2 log 2 p 2 = -A + log 2 A + - A_log 2 A_. (9) 

Since C has the same range for its values and is monotonically related to E(^), 
the concurrence can be regarded as a measure of entanglement in its own right. 

Employing local £7(2) x £7(2) rotations our \^f) can be transformed to the (4) 
Schmidt form i.e. we have 

l*)= E V^7l«>i^«)|/5> 2 %, tfU = VW = I, (10) 

j,a,(3=0,l 

where Ao,i = A +) _. Using this expression one can check that (see also [9]) 

-1-C = UDV T , Pl = UD 2 U\ p 2 = VD 2 V\ (11) 
v2 

where D is the diagonal matrix containing the square root of the eigenvalues X± in 
its diagonal. Taking the magnitude of the determinant of the expression for C in 
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Eq. (11) shows that local transformations presereve the concurrence, hence they do 
not change the degree of entanglement. This observation fulfills our expectations 
that entanglement can be changed only by global (i.e. U(4) transformations). 



III. Entanglement and the quaternionic Hopf fibration 

According to Equation (5) the set of normalized states is characterized by the 
constraint |a| 2 + \b\ 2 + |c| 2 + \d\ 2 = 2 (recall our convention of pulling out a factor 
of from the expansion coefficients of is the seven-sphere S 7 . The basic 

observation of [12] is that for understanding the geometry of two-qubit entangle- 
ment it is useful to fibre S 7 over the four dimensional sphere S 4 by employing 
the second Hopf-fibration. Moreover, it is convenient to introduce quaternionic 
notation for our basic quantities since the geometry of this fibration then easily 
described. Let us represent an element of S 7 by the quaternionic spinor, i.e. let 



1 

71 



b 

c 



U0 )=^( a + b j). (12) 



The quaternionic units i, j and k with squares equal to —1 satisfy the usual 
relations ij = — ji = k plus similar ones obtained by employing cyclic permutations 
of the symbols ijk. In this way an arbitrary quaternion q e H can be expressed 
as a pair of complex numbers 

q = qi + <?2i + <? 3 j + <?4k = (gi + q 2 i) + (<?3 + <?4i)j, (13) 

where the components qi, I = 1,2,3,4 are real numbers. The conjugate quater- 
nion q is obtained by changing the signs in front of the terms containing i, j and 
k 

q = qo - 9i i - <?2j - 93k = (q - 9i i) - (92 + ?3i)j- ( 14 ) 
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Recall also that due to the noncommutativity of quaternionic multiplication we 
have qp = (p)(q), and the norm squared of a quaternion q is defined as the real 
number \q\ 2 = qq. Now we are ready to define the second Hopf fibration by the 
map 7r as follows 

n:S 7 ^S\ (^j h- (2 Wl ^,K| 2 -M 2 ) = (£,£o), (15) 

where 

£ = (£1 + £ 2 i) + (£s + £ m £m = £o + ££ = 1, f* = 0,l,...4 (16) 

are Cartesian coordinates for .S 4 (summation for repeated indices is understood). 
Since 2u\Uq = (c+dj)(a + 6j) = (c+dj)(a— j&) — z + wj the Cartesian coordinates 
for S 4 can be expressed in terms of the quantities defined in (6) characterizing two- 
qubit entanglement, i.e. we have 

£ Q = ±y/l-\z\ 2 -\w\ 2 , Z 1+ & = z, &+& = w. (17) 

The basic result of [12] was that the mapping rr is entanglement sensitive. In 
this formalism this result is easily reproduced by noticing that submanifolds of 
fixed entanglement are characterized by fixed values for the concurrence C = 
Vl&iP + l^l 2 - H ence separable states are mapped to points of S 4 with vanish- 
ing values for the coordinates £3 and £4, i.e they are on the two-sphere £ 2 C S 4 
described by the constraint £q + £1 + £2 = 1- F° r maximally entangled states we 
have C = 1, then for these states we have £0 = £1 = £2 = 0. These states are 
parametrized by a great circle of the "equator" (which is a three-sphere S 3 ) of >S 4 . 

It is clear from Eq. (15) that multiplication of the quaternionic spinor from 
the right by a unit quaternion q (i.e. a quaternion with unit norm \q\ 2 = qq = l) 
leaves the coordinates ^ invariant. Since quaternions of unit norm ("quaternionic 
phases") form the group Sp(l) ~ SU(2) it means that entangled states related by 
an SU(2) rotation project to states of the same concurrence. This gauge degree of 
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freedom corresponds to the fiber of the second Hopf fibration. Since SU(2) ~ 5 3 
the fibration 71 is a one with total space 5 7 , base space 5 4 , and fiber 5 3 . The 
important result of the present paper is that the local gauge transformations of 
the Hopf fibration associated with this Sp(l) ~ 577(2) degree of freedom give rise 
to a geometric interpretation of local transformations in the second subsystem not 
changing the entanglement properties of our two-qubit system. The information 
available for the observer of the first subsystem is parametrized by the base space 
5 4 of the fibration as can be seen from the (7) form of the reduced density matrix 
p\. By exchanging the parameters b and c in the (12) definition we obtain another 
representation of the Hopf fibration with Cartesian coordinates 77^ for the corre- 
sponding four-sphere 5 4 . The assignment in this case reads 77 = ( + w j with ( 
defined by Eq. (6). In this case as was explained in Ref. [12] the roles of the two 
qubits are exchanged. Now the local gauge degree of freedom associated with our 
ignorance of the details of the first subsystem is represented by the fiber degree of 
freedom. The base space parametrizes the reduced density matrix p 2 . 

The relationship between base and fiber (i.e. the twisting of the bundle) is 
just the entanglement of the two qubits. A natural way of describing this twisting 
is via the means of introducing a connection on our bundle. Luckily for the second 
Hopf fibration we have a canonical connection the properties of which has been 
described in many places (see e.g. [17,18]). This connection is equivalent to the 
instanton connection well-known to physicists. Moreover, it can be related to a 
metric on HP 1 ~ 5 4 which is the quaternionic counterpart of the complex Fubini- 
Study metric on CP 1 ~ 5 2 . Our next task is to describe these quantities, and 
relate them to our basic ones of Section II. describing the phenomenon of two-qubit 
entanglement. 



IV. Sections, connections, and metrics 

First we introduce for two quaternionic spinors \v) and \u) the scalar product 
{v\u) = v a u a = vqUo + v\U\. Notice that right multiplication of our spinors with 
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the nonzero quaternion q yields the expression (vq\uq) = q{v\u)q. The vector space 
H 2 of quaternionic spinors with this scalar product is an example of a quaternionic 
Hilbert space (see Ref. 18. and references therein) . States in quaternionic quantum 
mechanics based on the space are represented by points of the quaternionic 
space of rays which is just HP^" 1 the N — 1 dimensional quaternionic projective 
space. It is amusing to see that the N = 2 case of interest for us yields the space 
HP 1 ~ S 4 , which is the quaternionic analogue of the usual Bloch-sphere represen- 
tation of complex spinors, i.e. we have CP 1 ~ S 2 . The complex Bloch-sphere is of 
basic importance for the geometric description of complex superposition, likewise 
the "quaternionic Bloch-sphere" plays a similar role for the geometrical description 
of quantum entanglement. Though this correspondence between entanglement and 
quaternionic quantum mechanics is an interesting idea to follow in its own right, 
here we work out merely the simplest N = 2 case and regard the quaternionic 
Hilbert space formalism merely as a comfortable representation. 

What is interesting for us is that two important geometric quantities can 
be defined on the space of normalized quaternionic spinors S 7 which pull back 
naturally to the base space HP 1 ~ S 4 , the metric and the connection. The first 
of these is related to the transition probability and the second to the 

transition amplitude (v\u) [18]. 

Indeed an invariant distance 



between two not identical, nonorthogonal quaternionic spinors representing en- 
tangled states can be defined. It is related to the distance along the geodesic 
connecting the points tt(\v)) and tt(\u)), representing the corresponding states in 
HP 1 , with respect to the metric which is the obvious quaternionic generalization 
of the well-known Fubini-Study metric. By using local coordinates Xk, k = 1, 2, 3, 4 
on an open set U C HP 1 parametrizing our spinor \u) and putting dl = A it is 
defined by the relation 




< A < 7T 



(18) 
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dl 2 = g k mdx k dx m = 4 [1 - \ (u + du\u) | 2 j . (19) 

Here it is understood that this equality is valid only up to terms higher than second 
order in the change of local coordinates. 

Moreover, since our entangled states are defined up to right multiplication 
with a unit quaternion it would be desirable to define a means for comparing 
the " quaternionic phases" of states of different entanglement. For the rule of 
comparing "phases" we adopt the definition that, two such states are "in phase" 
if (u + du\u) = 1 up to second order terms in du (the quaternionic analogue of the 
so called Pancharatnam connection [19] ). By introducing the quantity 



F = l-(u + du\u) (20) 

to be used later this rule can be restated as V = up to a second order term. 

In order to enable an explicit construction we have to chose a section for our 
bundle. This means that we have to adopt a choice for \u) G S 7 parametrized 
by points of S 4 . If the bundle is nontrivial the best we can do is to chose local 
sections. We chose the section 

with qq = 1. In this parametrization we have x = -ui(-uo) -1 hence it is valid on the 
coordinate patch U characterized by the constraint u ^ 0. Using this section we 
can pull back the (19) metric and (20) connection to the base space HP 1 yielding 
the formulae [18] 



2 _ Adxdx 

M ~ (i + M 2 ) 2 ' [ ] 



, Imxdx , 



n 



Here Imq = \{q — q) is the imaginary part of an arbitrary quaternion q. (Similarly 
the real part of q is defined by Keq = \{q + q).) The quantity 



A = Im 



(24) 



1 + Id 2 



is an sp(l) ~ s-u(2)-valued one-form (non-Abelian gauge-field) equivalent to the 
standard SU (2) instanton with self-dual curvature and second Chern-number = 
1 [17-18]. Notice that according to Eq. (20), Pancharatnam connection (r = 0) 
yields a condition for parallel translation of quaternionic phases. Indeed, using Eq. 
(23) with a suitable boundary condition we obtain the usual differential equation 
of parallel transport. For a curve C lying entirely in U with initial and end points 
being q(0) = 1 and q(r), we obtain the standard path ordered solution 



V Jc J 

Observe that our coordinates x = x\ + £21 + 2^+2^ used in the (21) section 
are related to the Cartesian coordinates £ = £i + £21 + £2j + and £0 as 



1 - 


X 


2 


1 + 1 


x\ 


2 



(26) 



Indeed, the coordinates (xi, X2, X3, X4) are obtained from stereographically pro- 
jecting the sphere S 4 from its south pole to R 4 U {00}. It is straightforward to 
check in these coordinates that dl 2 = d^ 2 + <i£ 2 + <i£f + <i£| + <i£| (the standard 
line element on S 4 ) is just (22). Let us define R 2 = £q + £ 2 + £|, then we have 
the relation 0<C 2 = 1 — R 2 < 1 where C is the concurrence. Using R and 
the polar coordinates (R, 6, $) on the unit ball B 3 originally parametrized by the 
coordinates (£o,£i,£2) we have 



d£o + ^1 + £3 = dR2 + R2 dtt, where dfl = dG 2 + sin 6 2 rf$ 2 . (27) 




(25) 
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Hence in terms of the concurrence C (or alternatively R) and these polar coordi- 
nates we have for the line element on the base space S 4 the expression 



dl 2 = + C 2 d X 2 + (1 - C 2 )rfO = + (1 - R 2 )d X 2 + R 2 dtf, (28) 

where w = \w\e ta,rgw = Ce lx . For separable states we have C = 0, hence this line 
element is reduced to dQ the one for the two-sphere S 2 corresponding to one of 
our separate qubits in the base. For maximally entangled states we have R = 
then the line element characterizing our base qubit reduces to the one of a circle 
i.e. dl 2 = dx 2 ■ The other qubit in all cases is associated with the fiber of unit 
quaternions. The relationship between the two qubits associated with the base 
and fiber in all cases can be described by parallel transport with respect to the 
connection T . To gain some insight into this relationship as a first step we have 
to express the pull-back one form A in terms of our complex coordinates z and w. 

For this we combine Eqs. (17) and (26) to see that the quaternionic phase 
of x (i.e. the unit quaternion p = ) can be expressed as 

P = T7\ = /m 2 + m | 2 G S P^ ~ SU W- ( 29 ) 
\x\ ^\z\ 2 + \w\ 2 

Moreover, since the sum of the squared magnitudes of £ and £o equals one it is 
useful to represent them as 



2\. 


x\ 




1 + 


\x 


2 



1 - 


x\ 


2 


1+1 


x\ 


2 



sin6> = — pr, cos 6= — (30) 



with < 9 < n The parametrization £ = sin6>p, £ = cos6> with p e S s cor- 
responds to introducing polar coordinates for S 4 . However unlike for the usual 
parametrization we favour z and w more than 9 hence we express it in terms of 



these quantities as 9 = arcsin(A/|z| 2 + \w\ 2 ). In this parametrization containing 
quantities characterizing the entanglement properties of our qubits the section of 
Eq. (21) reads as 
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/cosf\ = j_( V^a/i-N 2 -!^! 2 \ 

Vsmfp/ 9 y/2 1 ,/l qz \ z \2 _ | w |2 »+aJ J g> KJ 

where ± (=f) corresponds to sections over the northern or southern hemispheres. 
Since p,q G Sp(l) i.e. they are quaternionic phases, the parametrization in terms 
of 9 and the pair (p, q) is of the same form as the well-known parametrization of 
a complex spinor associated with the Bloch-sphere. However, the second equality 
also shows the meaning of these parameters in terms of the entanglement param- 
eters. Comparing Eq. (31) with Eq. (12) we realize that on the open set U we 
can always chose a section for which our parameter b equals zero. For later use 
here we also remark that in this (b = 0) parametrization a formula between our 
complex parameters z, w and ( of Eq. (6) holds 



C=UlT y/l-\z\ 2 -\w\A ™ ,Z , = sm 2 e/2—^, r = w/z. (32) 
2 V / 1 + \w/z\ z 1 + \r\ 2 

Expressing our (24) instanton gauge-potential in terms of the complex coordinates 
z and w we obtain on U the expression 

1/. sw ~r / 7 \ Im( zdz + wdw + (zdw — wdz)\) 
A = -(1 - cos6)lm(pdp) = ^ v JJJ -. (33) 

2 2y/l±y/l - \Z\ 2 - \ W \ 2 

We note that for another coordinate patch V with u\ ^ we would obtain an 
Sp(l) gauge-transformed expression for A [18]. From Eq. (33) we see that for 
separable states (w = 0) A defined on the submanifold S 2 (the boundary of the 
unit ball B 3 ) of 5 14 has the form 



A = -(1 - cos0)d$, where $ = arg^, cos6 = ±a/1 - \z\ 2 (34) 

which is the U(l) gauge potential of a magnetic monopole with pole strength \. 
Hence we see that when moving entirely in the £ 2 ~ 9B 3 submanifold of separable 
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states the relationship between the qubit in the base and the one in the fiber is 
characterized merely by the possible occurrence of a U(l) anholonomy factor. For 
maximally entangled states we have |«;| = 1 i.e. w = e lx and \z\ — 0, hence in 
this case we have the one-form A = — |lm (^) = — \d\ living on the great circle 
of the equator S 3 of S 4 . States parametrized by the points of this circle belong 
to the Sp(l) ~ S s fiber. Parallel transporting an element q G Sp(l) along this 
circle with respect to this one- form A yields an anholonomy factor of —1 or +1 
depending on the winding number of traversals beeing even or odd. In this way 
we have obtained an alternative proof for the well-known fact that the manifold 
of maximally entangled states in CP 3 is Sp(l)/Zi2 — S 3 /7i2. (S 7 is also fibered 
over CP 3 with the U(l) ~ S 1 fiber corresponds now to our circle parametrized 
by the angle %.) For states of intermediate entanglement labelled by the values 
of 9 in the interval < 6 < ^ from Eq. (31) we have a mixing between the 
complex coordinates z and w. This will result in a more complicated pattern for 
the anholonomy properties, reflecting the richness of the entanglement possibilities 
for the qubits. An explicit example for this phenomenon will be given in Section 
VI. 

V. The geometrical meaning of the Schmidt decomposition 

According to Eqs. (22) and (23) a line element and the pull-back of a connec- 
tion can be induced on our space HP 1 ~ S 4 which can be sliced to submanifolds 
of fixed entanglement. Now we make use of these facts to give geometrical inter- 
pretation to the Schmidt decomposition for two qubits. In order to do this we have 
to characterize a special subclass of geodesies in S 7 that project to geodesies on 
S 4 . For this we consider a curve C = \u(s)) C S 7 . Using Eq. (19) for this curve 
we have dl 2 = 4 (l — \ {u(s + ds)\u(s))\ 2 ) up to terms second order in ds. Taylor 
expanding this expression a formula for dl 2 is obtained 



dl 2 = 4\\Q(s)u(s) || 2 ds 2 , Q(s)=I-P(s), P(s) = \u(s))(u(s)l (35) 
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where \u(s)) = -^\u(s)). Right multiplication by a unit quaternion q(s) leaves in- 
variant the projector Q(s), but \u(s)) transforms to \ii(s))q(s) + \u(s))q(s). How- 
ever, since Q(s)\u(s)) = we see that dl 2 is gauge invariant hence it can be 
used to define the length of the "shadow" curve n(C) in S 4 of an arbitrary curve 
C C S 7 . Moreover, notice that expression (35) is also reparametrization invariant. 
Now we can characterize geodesies in S 4 in the following way. Geodesies in S 4 
are those curves 7r(C) through n(\u(si))) and n(\u(s2)}) for which the following 
reparametrization and gauge invariant functional 

L[C] = 2 f * ds\\Q(s)u(s)\l CcS 7 (36) 

J Sl 

is stationary. The variation of a similar functional for the complex case and the 
derivation of the geodesies was already given in Ref. [20]. For the quaternionic 
case the same steps has to be taken with the important difference that quaternions 
do not commute so we have to be careful in grouping terms. However, since 
the variation SL[C] is a real number it can be represented as the integral of the 
real part of a quaternion depending on s. Luckily we can cyclically permute the 
quaternionic entries under the operation of taking the real part (it is just the 
operation of taking the trace when we interpret the quaternions as two-by-two 
matrices) so the derivation is a straightforward excercise of following the steps 
described on pages of 220-223 of Ref. [20]. The result is the following. Using 
gauge invariance we can find a solution \u(s)) with initial vector \ui) = \u(0)) to 
the geodesic equation on S r which is horizontal i.e. parallel transported along 
the geodesic tt(\u(s))) with initial point 7r(\ui)) in S 4 . Moreover, exploiting the 
reparametrization invariance it is affinely parametrized i.e. ||m(s)|| is constant 
along CcS 7 . Such affine parametrizations are unique up to linear inhomogeneous 
changes in s, for convenience we chose the parametrization for which ||tt(s))|| = |. 
In particular a geodesic C starting from \v) which is the horizontal lift of the 
shadow geodesic n(C) connecting ir(\v)) and some other point n(\u)) is of the 
form 
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s s 

\u(s)) = cos-|0i) + sin-|0 2 ), {</>i\</> j) = Sij, », J = 1,2 (37) 



where 



|0i> = b), |02>= ^O-cos^l^/sin^, (38) 



and 



|«0 = l«>iTTT7, |W| = cos^. (39) 

Notice that in Eq. (39) in accordance with our conventions the quaternionic phase 
multiply the state \u) from the right. It is now understood that as was claimed in 
Eq. (18) A is the geodesic distance between the points n(\v) and tt(\u). 

Having the geodesic distance at our disposal, let us now define the measure 
of entanglement as the distance between an arbitrary entangled state and the 
separable state nearest to it. This idea has already been proposed and illustrated 
in [11] for the case of two qubits using algebraic geometric methods on the state 
space CP 3 . This space can be regarded as the base for an abelian U(l) fibration 
of S r hence it gives rise to an alternative parametrization for submanifolds of 
fixed entanglement. However in contrast to [11], when using instead the Hopf 
fibration of S 7 the fiber is the non-Abelian group Sp(l) of quaternionic phases 
making it possible for the two qubits to reside in different spaces, the base and the 
fiber respectively. The relationship between the qubits, i.e. their entanglement is 
measured by the twisting of the bundle. Hence it is instructive to see by giving an 
alternative proof, how naturally this measure of entanglement is encoded into the 
structure of the Hopf bundle. 

In order to see this first we chose a representative \u) of our entangled state 
in the (31) form with q = 1 and for the unknown separable state in the similar 
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form. This representative has b = d = hence C = 0, moreover it is already of the 
form of our standard section valid on U. We have to find the nearest separable state 
to | it), meaning we have to determine a and ip as a function of the entanglement 
coordinates z and w. In order to do this we have to maximize the expression 
cos 2 A vu /2 = \ (v\u}\ 2 . A short calculation yields for this quantity the expression 



| (v\u) | 2 = cos 2 — = cos 2 — cos 2 - + sin 2 — sin 2 - H — sinaRe(e t(p (z + wj)). (41) 
2 2 2 2 2 2 

Since our sections are merely local ones living on U we should exclude the south 
pole (<t = 9 = 7r) hence we have < 6 < n and < a < it. For these values (41) is 
maximal if (p = axgz . In this case we are left with the expression cos 2 |-cos 2 | + 
sin 2 1- sin 2 1 + ^\z\ sin a to be maximized with respect to changes in a. As one can 

I z I 

check this quantity is maximal provided tana = Hence separable states \v) 

nearest to our entangled state \u) labelled by the complex numbers z and w are 
characterized by the angles 





z 


2 


1 - 


W 


2 



coscr = ±Wl — (p = axgz. (42) 



Using these angles in Eq. (41) we obtain for the distance A uv the important 
formula 



cos 2 ^ = ^(l + v / T^). (43) 

Comparing our result with Eq. (8) we see that the value appearing in (43) is 
precisely the eigenvalue A + of the reduced density matrix also appearing in the 
Schmidt decomposition. Moreover, it is easy to see that the distance of our \u) 
from the state \v') orthogonal to \v) (this state is antipodal to \v) in S 4 ) is related 
to the other eigenvalue A_. Hence the (9) von-Neumann entropy as a measure 
of entanglement is just a special combination of lengths for the shorter and the 
longer segments of the geodesic linking our entangled state rc(\u)) to the surface 
>S 2 C S 4 of separable states. 
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What is the meaning of the separable state \v) nearest to the entangled state 
| it) ? It is just the quaternionic representative of one of the states \j) 1 (j = 0, 1) 
in the (10) Schmidt decomposition. In order to see this we have to diagonalize the 
reduced density matrices p\ and p 2 of Eq. (7). First we write these matrices in 
the form 



where the coordinates £ M are defined in Eq. (17), similarly 77^ is defined by the 
other Hopf fibration with the roles of b and c in (12) exchanged. Note also that 
from the five components of these vectors only the first three is used. We introduce 
the notation for these vectors 

/ Co \ / cos a \ ( 770 \ / cos t \ 

v = I Ci J = ]v| I sinacosi/? , t = I 771 J = |t| I sin r cose , (45) 
\ £2 / \ sin a sin <p J \ 772 / \ sin r sin e J 

where a and <p turn out to be precisely the quantities of (42), and the ones r 
and e are defined by replacing in (42) the coordinate z by ( of Eq. (6). In order 
to see this notice that v and t are elements of the unit ball B 3 i.e. we have 
|v| = |t| = a/1 — \w\ 2 < 1, and £ = cos6> = ±a/1 — \z\ 2 — \w\ 2 . In order to 
diagonalize the reduced density matrices we have to diagonalize vcr and ta. It 
is well-known that there is no global only local diagonalization of this problem 
over B 3 . The reason for this is the fact that the eigenstates of these operators 
form a nontrivial fibration over B 3 related to the first (complex) Hopf fibration. 
Local sections are well-known from studies concerning the geometric phase hence 
we merely refer to the result [19]. Eigensections of p\ and p 2 corresponding to the 
eigenvalue A + of (8) that are singular on the — £0 and — 770 axis are of the form 

w = U°P0 |t> = te^)- < 46 > 
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For the eigensections belonging to the eigenvalue A_ we have 

Looking at Eq. (46) we immediately see that the first vector |v) regarded as a 
quaternion is precisely \v) nearest to \u). 

Now the question arises: what is the meaning of |t) the Schmidt vector 
representing one of the orthonormal vectors \j) 2 (j = 0, 1) used for the other 
qubit? In the quaternionic notation we know that this vector resides in the fiber 
of the Hopf fibration, hence we shoud be able to recover it from the holonomy 
of our connection. In the following we will show that the representative of |t) 
is the unit quaternion Q obtained by parallel transporting our entangled vector 
\u) E S 7 with respect to the instanton connection along the geodesic segment 
between ir(\u)) and n(\v)) representing the closest separable state. 

In order to prove this first we represent our entangled state in the (31) form 
with q = 1, hence its matrix ^= and the diagonal matrix D of Eq. (11) has the 
following form 



C ( cos f \ /cos# 



f= . smf . * sinf —g D =[ o 2 sin^j- (48) 

Likewise using (46) and (47) the unitary matrices U and V diagonalizing on U the 
reduced density matrices p\ and p2 are of the form 



/ cosf -sinfe / cosf -sinfe^N 

V sin f e ^ cos f /' ~\sin^e ie cosf )' K ' 

(Recall that <p = axgz and e = arg£, Eq. (42) and a similar one with z replaced 
by C-) Putting these matrices in the first of Eq. (11) we obtain the relations 



-in T -e- = SinfS r f ,} W ' A , (50) 



2 cos y a/1 + \w/z\ 
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and 



r cos f cos 7T sin f sin ^ 1 

cos - = + — Sr A ^= • ( 51 ) 

2 cosf- cosf^ a/I + IW- 2 ! 2 

to be used later. A quick check for the phase of (50) is given by relation (32) 
giving e = arg ( = arg w/z which is confirmed by Eq. (50) too. 

Let |$i) = \v) where \v) is the (40) separable state a representative of the 
states over n(\v})\ We know that the point ir(\v)) is the nearest one to n(\u)). 
Consider now the (37) unique horizontal geodesic passing through the antipodal 
separable states \v) and \v') where the state |$ 2 ) = \v') is defined by Eqs. (38-39). 
This horizontal geodesic in S 7 is also passing through the state \u) which 
is apart from a quaternionic phase is our entangled state we have started with. 
Explicitly we have the relation 



. , (u\v) A, . . A. ,. . 

l " > niiT= cos 2 |, ' >+sin 2 |c> - (52) 



Multiplying this equation from the right by the quaternionic phase 



(v 


u) 


\(v 


u)\ 



and 



noticing that cos ^ and sin A are just the square roots of the eigenvalues of the 
reduced density matrices we obtain the form which looks like the quaternionic 
version of the Schmidt decomposition 



|«> = v^l^Q+v^Tl^P, Q=Prr r (53) 

\{v\u}\ 

where P is the quaternionic phase transforming \v') to \v ) having the (47) stan- 
dard form. In order to show that it is indeed the Schmidt decomposition we have 
to show that the quaternionic phase Q is somehow representing the other qubit 
in the Schmidt decomposition belonging to the fiber. In other words we have to 
show that the vector |t) of (46) corresponds to Q. Using the (53) definition for Q, 
the sections (31) and (40) with q = 1 and relations (50-51) one readily obtains 
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^ cos % cos % sin f sin § 1 + (w/z)\ r . r • . „, 

Q = + u \ / o = cos - + sin -e l£ j. (54) 



cos 



cosf y/l + \w/z\ 2 " 2 2 



Since according to Ref. [12] a generic state \i^q) in the fiber is defined as Q = 
IV'q) = c o|0)q + ci|1)q where the orthogonal states |0)q and |1)q are related to 
the choice Q = 1 and Q = j respectively, the correspondence between Q and |t) 
is established. Alternatively the reader can check by calculating \v)Q and then 
identifying the parameters a, 6, c, <i using (12) and (5) that this quaternionic spinor 
represents the separable state 



(cos ||0)i + sin le^ll)!) ® (cos^|0) 2 + sin^e ie |l) 2 ) = |v) ® |t). (55) 

Based on these considerations we have the following geometrical representa- 
tion of the Schmidt decomposition for two qubits. Chose first a special entangled 
state represented by a quaternionic spinor in the standard form (i.e. by putting 
q = 1 in Eq. (31)) and call this \u). This state is represented by a point u = 7r(|w)) 
in HP 1 ~ 5" 4 . Find the point v in the separable sub manifold £ 2 C 5" 4 nearest to 
u connected by the unique geodesic segment. (If \u) is not maximally entangled 
then we have a unique solution.) As a next step parallel transport \u) along this 
geodesic with respect to the instanton connection to the fiber over v. Then one 
pair from the biorthogonal Schmidt states is recovered as the standard (i.e. of 
the (40) form) section \v) over v and as the quaternionic phase Q between this 
section and the state obtained by parallel transport. The real expansion coefficient 
multiplying this pair forming the separable state \v)Q is just cos ^ where A is the 
geodesic distance between u and v. Repeat the same process for the antipodal 
point v 1 - of v in S A to obtain the the other pair of Schmidt states with expansion 
coefficient beeing the corresponding geodesic distance. 

If we use the Hopf fibration with the first qubit belonging to the base and 
the second to the fiber the gauge degree of freedom is manifested in the freedom 
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for the choice of the unitary matrix V in Eq. (11). This corresponds to a different 
choice for the local base belonging to the second qubit. In our representation 
for an arbitrary entangled state it means that we have a quaternionic spinor \u) 
this time with some fixed q ^ 1. Of course this new choice will not change our 
geometric interpretation of the Schmidt decomposition. It is easy to see that 
the new state representing the first qubit in the Schmidt decomposition is the 
same, and the second is obtained from Q by right multiplication by q. Hence we 
have proved that the Schmidt decomposition for two qubits is amenable for a nice 
geometric interpretation in terms of the anholonomy of the canonical connection 
on the Hopf fibration. Instantons were originally introduced as classical solutions 
of SU(2) gauge theories in Wick rotated space-time. They also play an important 
role in quantum field theories describing tunnelling between different vacua. It is 
amusing to find them here as the basic entities describing a fundamental aspect of 
quantum theory, two qubit entanglement. 



VI. Geometric phases 



We have seen that the instanton connection on the Hopf fibration plays a 
vital role for the geometrical description of entanglement for two qubits. A further 
interesting possibility to explore is to look at the non-Abelian (an)holonomy of the 
instanton connection, and reinterpret the results in the language of entanglement. 
For this purpose we have to somehow generate closed curves in the base manifold 
S 4 and calculate the quaternionic phases picked up by an initial quaternionic 
spinor after completing a circuit. Reinterpreting these spinors as entangled states 
the result of this non-Abelian parallel transport will be some final entangled state 
with very different form but the same concurrence. For different loops we obtain 
different anholonomy matrices, that can serve as quantum gates. This process is 
called anholonomic quantum computation [14]. 

However, S 4 is sliced to submanifolds of fixed entanglement, and we know 
that for separable states the anholonomy of the connection is Abelian. Therefore 
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interesting curves generating non-abelian anholonomy are the ones not restricted 
to the S 2 submanifold with concurrence C = 0. Moreover, a look at the (33) 
expression of A we see that the states characterized by the condition \z\ = 
describe another two-sphere S 2 C >S 4 , with another monopole-like gauge- field on 
it. According to Ref. [12] states with \z\ = are the ones with trivial Schmidt 
decomposition. For curves lying in the submanifold of these states the holonomy 
is again Abelian. 

Let us suppose then that we have a curve C lying entirely in W C S 4 not 
belonging to any of the submanifolds described above. In this case we have to 
calculate the quaternionic phase q[C] as given by formula 



q[C] =-p e ~$c A . (56) 

A more managable form for q[C] can be given by dividing the loop C into N 
segments characterized by the N points £o, £i> • • • Cat = Co £ U C S 4 . These points 
represent a whole family of gauge equivalent entangled states belonging to the fiber. 
Let P n = P(^n) = \ u (£n)) ( u (Cn) \ be the projectors represented by the quaternionic 
spinors \u n ) parametrized by the coordinates £ n . It is clear that P n is independent 
of the choice of representatives. By using (u(£ n+ i)\u(£ n )) ~ / — A(^ n ) + . . . one 
can show that [21] 



K£o)>Pe 



-§c A = 



hm P(£ Ar )P(£ iV _ 1 )...P(£ 1 )K£ C) )}. 



(57) 



Notice that Eq. (57) can also be understood in the following way. Given a curve 
C and its division by N points we can approximate it by a geodesic polygon. If we 
have an initial vector \uq) over the point £o we can parallel translate this vector to 



the fiber over the next point £i obtaining the vector \u\) 



{u 1 






«o)| 



. After the next 



(u 2 


Ul}{U! 


1*0 ) 


\(U2 


u 1 )(u 1 





limit iV — > oo by virtue of K^n+il^n)! = 1 — \dl 2 + ... we get the path ordered 
product of projectors verifying Eq. (57). 
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Hence in order to describe curves in S 4 we have to characterize a one pa- 
rameter family of rank one quaternionic projectors P(t). These are two-by-two 
quaternion-Hermitian matrices consisting of a single dyadic product of quater- 
nionic spinors. It is easy to see that such projectors are of the form 

P(t) = I(/ + r M £ M (t)), (58) 

where ^ G S 4 and the 2x2 quaternion-Hermitian matrices T^^ji = 0, . . .4 are 
defined as 



In order to prove this note that {T M , IV} = 215^ , hence P 2 = P. One can also 
prove that we have 



(u\T^\u)=^ where \u) = 




which is just another way of describing the (15) quaternionic Hopf fibration. Using 
this we have (u\P\u) = 1, i.e. P = \u)(u\. Using the (21) section one can obtain an 
explicit formula for P in terms of x. By virtue of the (26) correspondence between 
£ and x we get Eq. (58) as we have claimed. 

We know that the anholonomy of the instanton connection representing the 
parallel transport of entangled states is described by Eq. (57) with given 
by (58). However, it is desirable to simplify (57) by restricting our attention to 
a subclass of curves of the form P(t) = U(t)PW(t), where P is a fixed projector 
representing the initial entangled state in 5" 4 . U(t) is a one-parameter family of 
2x2 quaternion unitary matrices, i.e. U(t) G Sp(2) ~ Spin(b) (Spin(5)/Z 2 ~ 
50(5)). Such matrices are generated by quaternion skew-Hermitian matrices S 
in the form U(t) = e tS , where S can be expressed in terms of the generators of 
Spin(5) as S = ot^S^ with = |[r M , IV] and a^ u = —a V[li are real parameters. 
For a closed loop we should have P(2n) = P(0). 
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Now suppose we have an initial entangled state \u) £ S 7 corresponding to 
the projector P = \u)(u\ and a closed loop P(t) = U(t)PU^(t) satisfying P(2tt) = 
-P(O) = P. This loop is defined by some choice for the real parameters a^ v . Then 
we have the formula ( Proposition 7.4 of [22]) 

Jim PfoOPfor.O . . .P(£i)P = e ts (cos(t\\PSP\\)P - ggj^g^ Pflp) 

(61) 

where P = P(£o) 5 an d the quaternionic matrix norm is \ \B\\ = Trn(B'B). Using 
this result we obtain our final formula 



\u)Pe-*o A = e * s (cos(t||PSP||) - ^f^^) l«> 



(62) 



where \u) = |w(£o)) is the quaternionic spinor representing our initial entangled 
state. 

As an explicit example we take 



S = ifi (cos kT 3 - sin kT 4 ) = ^a M „ [r M , T v ] 0<k<2tv 



(63) 



meaning the real parameters are chosen as a\z = cosk and CK14 = — sin^. One 
can check that AS 2 = -I hence U(t) = e ts = cos(t/2) + sin(*/2)2£. Using this 
we have U(2ir) = —I meaning P(2-7r) = P(0) for any initial projector. Hence for 
changing values for k we obtain a parametrized family of loops C K . 

Let us chose the initial state to be the maximally entangled Bell-state 



i (i\ 



|M> = vl I j ) " w = 75 (|0>I |G>2 + 11)1 8 11>2> ' (64> 

Since this vector is the eigenvector of the matrix 1?3 the initial projector is P = 
P(£o) = ^(1 + ^3,). This means that the initial vector £ e S 4 for our loop C K has 
components = (0, 0, 0, 1, 0). Straightforward calculation shows that 
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PSP = -\smK[ ] f 1 ), \\PSP\\ = JsinK. (65) 



4 V i -k J ' 11 1 2 
Putting these results into Eq. (62) we get the result 



1 Z' 1 A p ~§ C A _ ( cos(7rsinK) i sin(-7r sin k) \ 1 fl 
y/2 \j J 6 ° ~ V i sin(7r sin k) cos(tt sin «) J ^ \ j 

Hence after parallel transporting the Bell state \u) along C it will pick up a quater- 
nionic phase of the form 

q[C] = Pe $c A = — cos(7r sin k) — sm(n sin «;)k. (67) 

From this equation and the (12) definition one can read off the complex coeffi- 
cients identifying the new entangled state. These are a = d = — cos(7rsinft) and 
b = c = — i sin(7r sin k). Since C — \ad — bc\ = 1 we see that the anholonomy for our 
families of closed loops C K is not changing the degree of entanglement. This is not 
surprising since as the reader can check the SU(2) ~ Sp(l) anholonomy transfor- 
mations associated with an arbitrary loop C belong to the local transformations 
manipulating only the second qubit. 

For the choice k = ^ we have q[C] = — k. In this case the new state is 




\u)q[C] = ( • ) ^ W) = (|0>i ® |1>2 + ® |0> 2 ) (68) 

hence for the loop C^/q we obtain up to a sign another Bell state. 

We have not clarified the nature of our loops C K , < k < 2tt yet. In order 
to do this we write out explicitly U (t) as 



U(t) 



cos I + sin | (j cos k — k sin k) 






cos I — sin | (j cos k — k sin k) 



(69) 



and calculate the matrix 
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Straightforward calculation shows that the family of curves C K in S 4 is given by 



/ \ 

sin t cos k 


cos 2 1 — sin 2 1 cos 2k 



(71) 



\ sin 2 1 sin 2k J 

Noticing that £ = we see that our family of loops lies in the equator of S 4 which 
is a three-sphere S 3 . Moreover, recalling Eq. (17) we see that z(t; k) = sin t cos k 
and w(t; k) = cos 2 | — sin 2 |e~ 2m . Using C = \w\ the concurrence as a function of 
t and k can be determined 



C = (l + cos 2 t-sin 2 tcos2K). (72) 

We see that the loops visit sub manifolds of different entanglement, except for 
the case k = ^, z = when the loop degenerates to the starting point. It can 
be understood as follows. Our U(t) is an element of the subgroup Spin(4) C 
Spin(5) ~ Sp(2). Spin{4)/7i2 — SO (A) is the four dimensional rotation group. 
The terms Til^ and r^r^ in (63) generate 50(4) rotations in the 13 and 14 
planes respectively. Clearly for k = ^ we have merely 14 rotations not changing 
the £3 = w = 1 constraint characterizing our initial Bell state. On the other hand 
for k = we have 13 rotations. In this case we have z(t) = sint and w(t) — cost 
hence C = | cos£|. This loop starting from the maximally entangled submanifold 
then crosses the separable surface at t = \ meeting the maximally entangled 
surface again at t = n and gets back to the separable surface at t = ^ and then 
to the initial state at t = 2n. The anholonomy matrix for this curve from (66) is 
—I reproducing the sign change of quaternionic spinors under a 2tt rotation. For 
the case « = f studied in Eq. (68) we get for the concurrence C = ^V3 + cos 2 t 
along the loop. 
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Since 5 4 ~ SO(5)/SO(4) (i.e. 50(5) acts transitively on 5 4 ) we can also 
consider more general loops generated by SO (5) rotations. Such loops are gen- 
erated by an S also containing terms of the form aoiToTk^k = 1,2,3,4. There 
is also the possibility of considering loops generated by ordinary SO (3) rotations. 
However, the structure of 50(3) orbits on 5 4 is more complicated. For the details 
of these orbits see Ref. [23], here we merely give the (skew-Hermitian) 50(3) 
generators needed to generate such loops 



One can check that the usual relations [Jj,Jk] = tjkmJm hold. It is now clear 
that using an explicit form for S and employing formula (61) one can calculate 
the anholonomy of an arbitrary entangled state represented by a fixed projector 
P for loops regarded as SO(n) orbits (n = 3,4, 5) of the form P(t) = U(t)PW(t). 

It is now obvious that if we have a means of generating a prescribed set 
of closed loops on HP 1 ~ 5 4 through an entangled state of fixed concurrence 
then we can associate to this loop an SU(2) anholonomy matrix. Moreover, this 
correspondence between loops and anholonomy matrices can be useful in building 
one-qubit gates in the spirit of anholonomic quantum computation [14]. Due to 
the geometric nature of these gates, quantum information processing is expected 
to be fault tolerant. 

With this possibility in sight the question arises: how to generate loops 
on HP 1 ? The first possibility is the case of adiabatic evolution. For this one 
can consider the parametrized family of 2 x 2 quaternion-Hermitian Hamiltonians 
H(^) = reinterpreted as 4 x 4 complex-Hermitian Hamiltonians with two 

Kramers degenerate doublets [24]. Note, however that unlike in Ref. [24] now 
we have to regard H(£) as a parametrized family of Hamiltonians coupling two 
qubits. In this picture we can view the entangled states as a parametrized family 
of eigenstates (eigensections) of these Hamiltonians. 
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In order to be more explicit we write H(£) in the form 



H(X) = X mn J m J n = ^ = 0,1,... 4, m, n= 1,2,3 (74) 

where 



/ 


-6 + 7^0 


-6 


"6 \ 




-6 


6 + 7^0 


£4 


V 


-6 


£4 





(75) 



In this way we have made a mapping from S 4 to the space of real traceless sym- 
metric 3x3 matrices satisfying |TrX 2 = 1 i.e. to the space of unit quadrupoles. 
H(X) are precisely the quadrupole Hamiltonians studied by Avron et.al in Ref. 
[24]. However, now these Hamiltonians have a different interpretation. We are 
not regarding these operators as describing a spin | particle in an adiabatically 
changing quadrupole electric field with the underlying Hilbert space beeing C 4 . 
We rather regard them as a parametrized set of coupling Hamiltonians for two 
qubits with Hilbert space C 2 x C 2 . This correspondence provides a nice formal- 
ism to label entangled states with the eigenstates of quadrupole Hamiltonians. 
For example separable states (£3 =£4 = 0) are represented with a block diago- 
nal quadrupole. Maximally entangled states are the ones with a complementary 
structure for X mn , i.e. £0 = £1 = £2 = 0. The Bell state of (64) is represented by 
the quadrupole with the only nonvanishing components X\z = X31 = —l/y/S. 

Now we fix a unit quadrupole representing an entangled state. Changing the 
quadrupole components (i.e. the coupling between the two subsystems ) then we 
adiabatically trace out a loop C in S A . Thanks to the adiabatic theorem [19] the 
time dependent Schrodinger equation maps an initial eigenstate into an instanta- 
neous eigenstate of the quadrupole Hamiltonian hence this state evolves along an 
open curve in S 7 (the set of normalized states in our Hilbert space H ~ C 2 x C 2 
) whose shadow in S 4 is C. The difference between the initial and final state in 
S 7 is an SU{2) matrix (an Sp(l) quaternionic phase) containing a dynamic and a 
geometric part. The geometric part in this case can be separated, it is precisely the 
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anholonomy matrix corresponding to the operation of our quantum gate. In this 
way we have implemented quantum computation via quantum adiabatic evolution 
in the slowly changing environment of unit quadrupoles. 

The other possibility for realizing closed curves on HP 1 is by the non- 
Abelian version of the Aharonov-Anandan phase [25] . In this case the assumption 
of adiabaticity is relaxed, by introducing cyclic solutions of the time dependent 
Schrodinger equation. These are such solutions of the Schrodinger equation in S 7 
whose shadow curves are closed in HP 1 . (In this case the solution curve is not an 
eigenstate of the instantaneous Hamiltonian.) In order to obtain an example of a 
cyclic evolution process in the quaternionic representation we have to find a 2 x 2 
quaternion unitary matrix U (t) and a T e R for which we have 

U(T)\u(0)) = \u(0))p, peSp(l). (76) 

We can specify such a U (t) by giving its quaternion skew-Hermitian generator S. 
A nice example of this kind is given by the choice 

S(t) = \[H(X(t)\H(X(t))} = \[H(r,(t)),H(f,(t))} = \\T„T v ] V ^t)v v (t), (77) 

for an n^(i) defining a closed loop (rj^T) = ?y M (0)) in another four sphere S^. 
Hence U(t) is a Sp(2) ~ Spin(5) rotation of the form U(t) = e s ^ with a^i) = 
(Vn(t)Vi>(t) — Vvit)^^))- It is important not to mix the four-sphere defining 
the parameters of the Spin(5) rotation with the other four-sphere HP 1 ~ 5" 4 the 
"space of states" for the entangled two qubits stratified into submanifolds of fixed 
entanglement. We can of course identify the parameters of with the space 
of unit quadrupoles (this is reflected in the first equality in (77)), but then the 
relationship between the quadrupole components and the coordinates of HP 1 is 
not canonical as was in the case of adiabatic evolution. 
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It is easy to show that the time evolution operator along the curve C in S 7 
is precisely the operator of parallel transport along n(C') = C C U C HP 1 i.e. we 
have [18] 



Pe -f c ,S(t)dt 



n — 1 2 1 

u(0)) = lim P(T)P( T)...P{-T)P{-T) |«(0)>, (78) 



where P(t) is the projector belonging to the quadrupole operator H(rj(t)) = 
H(X(t)) = T M ?y M (t) (the relationship between X mn and is of the same form 
as in (75)). Note however, that unlike in the adiabatic case |it(0)) is not an eigen- 
vector of the instantaneous Hamiltonian corresponding to S(t). The evolution is 
nonadiabatic and cyclic. Since the dynamical phase [19] is just the integral of 
(u(t)\S(t)\u(t)) = where \u(t)) are the instantaneous eigenstates of H(t), the an- 
holonomy of the evolution is purely geometric. Hence we see that quadrupole-like 
Hamiltonians are capable of generating closed curves via the standard Schrodinger 
type of evolution (both adiabatic and nonadiabatic) in the stratification manifold 
of two-qubits, enabling an implementation for anholonomic quantum computa- 
tion. The quantum gates obtained in this way are anholonomy transformations of 
the instanton connection our basic entity governing the entanglement properties 
of two qubits. 



We have seen that the base space of our fibration is S 4 which is parametrized 
by the reduced density matrix pi of Eq. (7). In this section we want to make this 
statement more precise. As we see from Eq. (44) the space of density matrices 
is the three dimensional unit ball B 3 . The relationship between B 3 and S 4 is 
clarified by Eq. (28) which shows that the line element on S 4 is a combination of 
a line element on B 3 of the form 



VII. Density matrices 



Adl 2 B = 



dC 2 



+ (1 -C 2 )rfO, 



(79) 



1-C 2 
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and the line element of a circle Sg, dl 2 sl = C 2 dx 2 with its radius parametrized by 
the value of the concurrence C. The line element in (79) is the one corresponding 
to the Bures metric [26] on the space of density matrices defined by 

d B (p, u) = ^2- 2Tv(u 1 / 2 pu 1 / 2 ) 1 / 2 , (80) 

where p and u> are density matrices. Indeed, by restricting our attention to 2 x 2 
density matrices for the infinitesimal form of (79) we obtain the relation [27] 

dl% = ]-Ti{dpdp) + d(Detp) 1/2 d(Detp) 1/2 . (81) 

Parametrizing p as in the first expressions of (44) and (45) and recalling that 
Detp = C 2 /4 we obtain (79). Hence the line element on 5" 4 is of the form 

dl 2 = Adl 2 B + C 2 d X 2 . (82) 

The space of 2 x 2 density matrices V 2 is stratified into the submanifold 
of rank one density matrices T> 2 (1) ~ <9B 3 ~ S 2 , and the submanifold of rank 
two density matrices T>2(2) ~ IntB 3 . We realize that rank one density matrices 
correspond to separable states with C = 0, and rank two density matrices to the 
entangled ones with C > 0. 

Having clarified the correspondence between reduced density matrices and 
the base space, we now turn to total space S 7 of the Hopf bundle. S 7 can be 
regarded as the space of purifications for the reduced density matrix p\ . The 
structure of these purifications was studied in [28]. Here we use these results to 
give some new insight to the geometry of two qubit entanglement. 

Let us denote by M.2 the space of complex 2x2 matrices with elements 
A G M.2 satisfying the constraint TrAA^ = 1. It is clear that in our case A = 
-^C of Eq. (1), moreover due to normalization M. 2 — S 7 . Then the mapping 
/ : A G M.2 !— > Pi = AA" 1 " G T>2 defines a stratification (see also Eq. (2)). This 
stratification is the union of fibre bundles 
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M 2 = \jM 2 (k)^V 2 = \Jv 2 (k), fc=l,2 (83) 
k k 

where M. 2 {k) denotes the manifold of rank k matrices in M. 2 . 

For the nonsingular case DetA = ^DetC = \w ^ 0. Write w = \w\e lx = 
Ce lx . Since the quantity AA^ is invariant with respect to right multiplication by 
an element of U{2) ( the unitary V of Eq. (11)) the mapping f 2 : M 2 (2) — > T> 2 (2) 
is a principal bundle with U{2) fiber. There is a subbundle B 2 of this bundle 
defined by the constraint x = 0, i.e. DetA e R+. As one can check B 2 is an 
577(2) ~ Sp(l) bundle over IntB 3 . This bundle will turn out to be important in 
the following considerations. 

For the singular case it was proved in Ref. [28] that we have Ai 2 (l) ~ 
5 3 X{/(i) 5 3 hence we have an imbedding of the complex Hopf bundle into M. 2 (l). 
This is just the precise mathematical statement of the observation in [12] that in 
this case the quaternionic Hopf fibration can be " iterated" to include the complex 
Hopf fibration. Moreover, in the light of this result it is not surprising that we 
have obtained in Section IV. the magnetic monopole connection (the canonical 
connection on the complex Hopf bundle) on this stratum. 

An interesting topic to discuss here is the relationship between Uhlmann's 
connection for parallel transport for mixed states [29] and our connection governing 
the entanglement properties of our qubits. Uhlmann introduced "amplitudes" of 
density matrices. In our case these are just the matrices Ai and A2 purifying two 
different nonsingular reduced density matrices p and u>. It is understood that p 
and uj correspond to the reduced density matrices of two different entangled states 
with coordinates 2:1,2 , ^1,2 and concurrences Ci j2 7^ 0. This means that we have 
p = AiA\ = ^C\C\ and u = A 2 A^ = \C 2 C\. According to Uhlmann two such 
amplitudes are parallel iff 



= C\C X > 0. 
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(84) 



It was also shown that on M.2(2) this parallel transport law is implemented by a 
connection form A satisfying 



CUC -dC ] C = C ] CA + AC ] C. 



(85) 



Following Ref. [28] it is straightforward to check that Uhlmann's connection A on 
the subbundle £> 2 is just our A instanton connection. First we take the trace of 
(85) then we use the identity MN + NM = Tr(M)N + MTr(N) + (Tr(MiV) - 
Tr(M)Tr(iV)) valid for our 2x2 complex matrices to get 



Here we have taken into account Tr(C^C) = 2 and the restriction Tr(^4) = 
which is valid on the subbundle £>2 [28]. Now writing out A = Im(u\du) as a pure 
imaginary quaternion-valued one form in the quaternionic notation of (12) and A 
as an s-u (2) -valued one form in the complex notation for C as given by (5), one 
can see that the two connections coincide. It is important to stress once again 
that A = A only for entangled states lying in the subbundle with x = argw = 0. 
According to (82) for such states the line element of the metric on S 4 is just 
the Bures line element so the identification of the relevant connections is not 
surprising. Since the connections for this class of entangled states coincide, so 
does the anholonomy properties of entangled states parallel transported in the 
subbundle £> 2 of M.2(2). Using the results of the previous section we can generate 
curves in the space of reduced density matrices and calculate the anholonomy of 
the entangled states regarded now as purifications lying in &2- In this way in 
the context of two-qubit entanglement we managed to find a nice application for 
Uhlmann's parallel transport for density matrices. 

In the following we would like to make some comments on the conformal 
structure of our metric on S 4 . Let us consider first the line element dl\ associated 
with the Bures metric. Use the notation R 2 = £ 2 + £i + £2 = 1 — C 2 and introduce 
the new coordinate —00 < f3 < 00 via the relation R = tanh/3. Since — 1 < R < 




(86) 



35 



1 this is a good parametrization for R. Notice that (3 behaves as the rapidity 
parameter in the special theory of relativity. For (3 = we obtain maximally 
entangled states and for (3 = ±oo we get the separable states. It is straightforward 
to check that using this new parametrization the Bures line element takes the 
following form 



sheet of the double sheeted hyperboloid H 3 we can conclude that the Bures metric 
is conformally equivalent to the standard metric of hyperbolic geometry. The 
hyperboloid H 3 can be stereographically projected to the Poincare ball B 3 with 
the standard Poincare metric on it. Since the space of 2 x 2 density matrices is 
just B 3 the Bures metric is up to a conformal factor is just the Poincare metric. 
Hence the space of density matrices can be given a hyperbolic structure. In this 
hyperbolic structure separable states are infinitely far away (i.e. on the boundary 
<9B 3 ) from the maximally entangled ones. This observation has already been made 
in a different context by Ungar in his study of gyrovector spaces and the geometry 
of 2 x 2 density matrices [30]. We can take one step more by realizing that the (28) 
line element for the four-sphere can also be given a conformal form after noticing 
that (1 - R 2 )dx 2 = dx 2 / cosh/3 i.e. we have 

dl 2 = \- (d(3 2 + sinh 2 pdfl + d X 2 ) ■ (88) 

cosh (3 

This means that conformally we have 



Notice that for this conformal equivalence we had to remove the boundary of B 3 
i.e. the separable states. We can also understand Eq. (89) by looking at the 
standard line element on R 4 , dl^ = dx\ + dx\ + dx 2 + dx 2 . The Xk 7 k = 1, 2, 3, 4 



Since the line element dl1, 3 = 



Adl 2 B = 




IntB 3 x S 1 . 



(89) 
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are just the stereographically projected coordinates of (26). Let r, and \ be polar 
coordinates on the X3 — X4 plane. According to (26) this plane is just the one 
corresponding to our complex coordinate w = £3 + i^4 = Ce %x . Then we can write 



( dv ~\~ dx ~\~ dx^ 1 

<^r 4 = r2 ( £ ~ + d X 2 J > < r < 00, -00 < xi,x 2 < 00. (90) 



The metric (dr 2 + dx\ + dx?,)/r 2 is the hyperbolic Poincare metric of the upper 
half space U 3 which can also be mapped to the B 3 Poincare ball [31]. The metric 
on U 3 is singular at r = (separable states), hence we have the identification 
R 4 — R 2 ~ U 3 x S . After stereographic projection the left hand side is conformally 
equivalent to the four sphere S 4 : and the right hand side by making use of the 
topological equivalence of U 3 and B 3 is just B 3 x S 1 . 

Closing this section as an interesting application of these ideas we calculate 
the geodesic distance between two entangled states and |<E») represented by 
the quaternionic spinors \u) and \v). These spinors map to points ^ and r\ v 
[i, v = 0, 1, . . . 4 of S 4 via the projection of the Hopf fibration. The first three 
components of the vectors £ M and r} v we denote by the vectors u and v. We relate 
as usual the remaning components to the concurrences as £3 + = C\e lXl and 
V3 + = C2e %X2 . We denote the reduced density matrices corresponding to \u) 
and \v) by g and uj. Then the transition probability related to the (18) geodesic 
distance between our entangled states is given by the formula 



I {u\v)\ 2 = Tr H (P u Pv) = ^(1 + ^) = \ (1 + uv + dC 2 cos( X i - X2)) , (91) 

as can be checked by using the (58) form of our projectors and the Clifford algebra 
properties of the V matrices. For the subbundle Ei characterized by the constraint 
Xi = X2 = this expression according to Ref. [28] is just the Bures fidelity 
[Tr^ 1 / 2 /^ 1 / 2 ) 1 / 2 ] 2 occurring in the (80) Bures distance. In terms of the rapidities 
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P u and f3 v after introducing the Lorentz factors 7 U)V = cosh/3 UjV = 1/ a/1 — |u, v| 2 
we can rewrite this expression for the distance of our entangled states as 

cos^ = \{u\v)\ 2 = — *— (7t +cos(xi - X2)), It = 7u7v(l + uv) (92) 
2 27 U 7 V 

where 7t is obtained by the addition law of velocities for u and v in special 
relativity (the cosine theorem of hyperbolic geometry [30]). As a special case for 
the subbundle B2 and for states with xi = Xi we obtain the formula of Ref. [32] for 
the Bures fidelity amenable for a nice interpretation in hyperbolic trigonometry. 

VIII. Conclusions 

In this paper we related the basic quantities of two-qubit entanglement to 
the geometrical structure of the quaternionic Hopf fibration. The entangled state 
was represented by an element of the bundle space S 7 . One of our qubits was as- 
sociated with the base (S ) and the other with the fiber (S 3 ) of this fibration. The 
nontriviality of the fibration i.e. the twisting of the bundle have been connected 
with the entanglement of the qubits via the use of the canonical connection and 
the natural metric on S 7 . These quantities pull back via the use of sections to 
the base giving rise to the instanton gauge field with self-dual curvature, and the 
Mannoury-Fubini-Study metric. 

The base space can be startified to submanifolds of fixed entanglement. Sep- 
arable states occupy a two sphere S 2 , maximally entangled states a great circle S 1 
of the equator of HP 1 ~ S 4 . The complement of separable states in the base can 
be conformally represented as JntB 3 x S 1 . Here IntB 3 is the space of rank two re- 
duced density matrices. The measure of entanglement was defined as the length of 
the shortest geodesic with respect to the Mannoury-Fubini-Study metric between 
the entangled state in question and the nearest separable state. This nearest sep- 
arabe state expressed in the form of the standard section is one of the Schmidt 
states appearing in the Schmidt decomposition of the entangled state. The other 
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Schmidt state is reproduced from the quaternionic phase between this separable 
state and the one obtained by parallel transport along the shortest geodesic from 
our initial entangled state. The other pair of Schmidt states is obtained by a sim- 
ilar procedure carried through by using the longer part of the geodesic, and the 
corresponding states. 

We examined the anholonomy properties realized by quaternionic geometric 
phases for the two-qubit entangled states. We have shown that for quadrupole-like 
Hamiltonians we can generate closed curves using the standard Schrodinger type 
of evolutions (both adiabatic and non-adiabatic) in the stratification manifold of 
two qubits. These evolutions enable an implementation for anholonomic quan- 
tum computation. For a specific family of curves we explicitly constructed the 
anholonomy transformations corresponding to a special class of quantum gates. 

By looking at the total space S 7 of our fibration as the space of purifica- 
tions for our reduced density matrices we have started working out a dictionary 
between the non-Abelian SU(2) geometric phase governing the entanglement of 
our qubits and Uhlmann's law of parallel transport for purifications over the man- 
ifold of density matrices. The space of such purifications forms a startification, 
i.e. a collection of fibre bundles. These bundles provide the natural geometrical 
setting for a deeper understanding of the results of Ref. [12]. We have shown how 
the Bures distance between reduced density matrices is related to the geodesic dis- 
tance between entangled states. We have also reformulated the known relationship 
between the instanton connection and Uhlmann's connection from the viewpoint 
of two-qubit entanglement. Based on results of Ungar et.al we made connection 
with the geometric data of the Hopf fibration, entanglement, and the hyperbolic 
geometry of the space of reduced density matrices. 

It would be interesting to see whether the usefulness of quaternions is merely 
a specific property of two qubit entanglement or it can be generalized for other 
entangled quantum systems. In the context of quaternionic quantum mechanics, 
or merely as a convenient representation we can consider the higher quaternionic 
Hopf fibrations n : S 4n+S — > HP" with again an Sp(l) fiber. The connection 
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and metric used here can be generalized and has already been discussed elsewhere 
[19]. In this case one can use n + 1 component quaternionic vectors somehow 
representing entangled states. However, the structure of these bundles should be 
only capable of incorporating the twisting of merely one qubit (belonging to the 
fiber) with respect to some other (possibly also entangled) state represented by the 
base degree of freedom. The other possibility to utilize the nontrivial fibre bundle 
structure as a geometric representation for quantum entanglement is to consider 
spin bundles , based on higher dimensional Clifford algebras. In particular the 
2 n dimensional spin bundles over 2n dimensional spheres seems to be promising. 
For the n = 2 case we get essentially the Hopf fibration. For the general case 
the fibre is SU(2 n ~ 1 ), and the tensor product structure of higher dimensional V 
matrices hints a possible use for understanding a special subclass of entangled 
systems in the language of higher dimensional monopoles as connections on such 
bundles. Apart from the fact whether these expectations are realized or not we 
hope that we have convinced the reader that the quaternionic Hopf fibration with 
its instanton connection is tailor made to unveil the basic structure of two-qubit 
entanglement. 
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